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Abstract. We demonstrate efficient and reversible mapping of a light field onto a 
thulium-doped crystal using an atomic frequency comb (AFC). Thanks to an accurate 
spectral preparation of the sample, we reach an efficiency of 9%. Our interpretation 
of the data is based on an original spectral analysis of the AFC. By independently 
measuring the absorption spectrum, we show that the efficiency is both limited by the 
available optical thickness and the preparation procedure at large absorption depth for 
a given bandwidth. The experiment is repeated with less than one photon per pulse 
and single photon counting detectors. We clearly observe that the AFC protocol is 
compatible with the noise level required for weak quantum field storage. 



PACS numbers: 42.50.Md, 42.50.Gy, 03.67.-a 
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1. Introduction 

The ultimate control of a quantum light field naturally demands a good understanding 
of its interaction with matter. In that context the reversible mapping of a field into 
an atomic system represents a fundamental challenge. The experimental quest of a 
quantum memory for light has produced many accomplishments using different systems. 
This domain has also been motivated by the realization of a quantum repeater for which 
the memory is a critical element pp. Benefitting from a long experience, most of the 
recent achievements have involved atomic vapors [21 El El E] . More recently atom- 
like impurities in solid state matrices have been considered as excellent candidates for a 
quantum interface. They can indeed have good coherence properties [HE]. One can take 
advantage of their complete absence of motion to demonstrate extremely long coherence 
time [HI Uni EH] • The rare-earth ion doped crystals (REIC) are particularly interesting 
because they are relatively easy to produce and commercially largely available. The 
most promising storage protocols fully use the primary feature of the REIC: the optical 
transition have a long coherence lifetime and a large inhomogeneous broadening. These 
schemes are in the lineage of the photon echo technique and are intimately related to 
a dipole rephasing. The direct excitation of the optical coherences over a wide spectral 
range (larger than the homogeneous width) is usually the initial storage step. The role 
of the protocol is to rephase the dipoles by externally controlling the detuning of the 
transition [T21 [131 HU 03] > which can lead to large storage efficiencies [T6] . 

More recently a protocol based on the preparation of an Atomic Frequency Comb 
(AFC) has allowed light storage at the single photon level [T71 [18]. The protocol is 
based on the absorption of a photon by a spectrally periodic structure of narrow peaks 
(a comb of periodicity 1/T). Right after absorption, each spectral class of dipoles starts 
dephasing with respect to each other. But at the time T, the accumulated phases are all 
equal to first order to an integer number of 2tt (0 modulo 2ir): this coherence build-up 
triggers a retrieval. Using a temporal image, this procedure can be interpreted as a 
dipole rephasing. Alternatively a spectral image can also be employed. The periodic 
spectral filter is a diffraction grating. It produces a delayed response in the time domain 
(the retrieval). This spectral image has been widely used to explain the stimulated 
photon echo (SPE). The AFC is in the direct lineage of the SPE where a periodic 
modulation of the absorption profile is considered. In that sense, the AFC represents 
an extreme situation: high narrow absorbing peaks are separated by fully transparent 
region. Such a comb can absorb close to 100% of the light and reduces intrinsic dephasing 
during the storage time 1/T. As a consequence the retrieval efficiency should be high in 
the forward direction (54%) and ideally perfect (100%) in the backward configuration 
[T7] . In order to achieve high efficiency, it is thus crucial to achieve a precise spectral 
shaping of the inhomogeneous absorption profile. The AFC is not only an interesting 
alternative storage protocol. It has been shown to have the highest temporal multimode 
capacity in principle [19] , which is a critical figure of merit in the prospect of quantum 
communication [20J. 
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In this paper we show how the AFC protocol can be efficiently implemented in 
an appropriate crystal namely Tm 3+ :YAG . Our work is focused on the initial step 
of the protocol where the signal is stored onto the optical coherences. Our results in 
terms of storage efficiency correspond to more than an order of magnitude improvement 
as compared to previous proof of principle realization performed in Nd 3+ :YV04 [18J. 
This improvement is possible because the Tm:YAG allows for a better preparation of 
the AFC due to an efficient optical pumping and a narrower homogeneous linewidth. 
More generally the material properties gives the ultimate width of the spectral selection 
(presence of spectral diffusion, superhyperfine effect ...). Additionally any external 
source of broadening as the laser linewidth can reduce the efficiency for a given 
bandwidth. This limiting factor is not present in our case since we use a stabilized 
laser. 

2. Efficient light storage in Tm 3+ :YAG 

Thulium has an interaction wavelength (793nm) easily accessible with laser diode as 
compared to other rare-earth ions (praseodynium and europium). This is particularly 
convenient when a complex spectral preparation sequence is required. The trivalent 
thulium is also known to have a long optical coherence time from 10 to 100 /is depending 
on the crystalline matrix, dopant concentration and magnetic field [2 lj . The spin 
coherence lifetime has been also observed to be long in the ground state [TT]. Very 
precise spectral tailoring will ultimately lead to a high number of temporal modes stored 
in the medium, which primarily scales as the number of spectral channels accessible in 
the absorption profile. The long hole-burning lifetime observed under an appropriate 
magnetic field will be especially useful for the initial state preparation [22]. This ensures 
an efficient population transfer between hyperfine levels (Zeeman states) by optical 
pumping [23] . 

2.1. Experimental setup 

To implement the AFC protocol in Tm 3+ YAG , we apply a 210G magnetic field along 
the [001] crystalline axis. This splits the ground and excited levels into a nuclear spin 
doublet by A g = 6MHZ and A e = 1.3 MHz respectively [24] (see fig. [Q. This is 
sufficient to create a few megahertz interaction bandwidth (~ 3 MHz in our case). Our 
0.5% doped- Tm 3+ YAG crystal is immersed in liquid helium at 2.3K. This permits 
a good coherence lifetime (T 2 = 30/xs) and then an accurate spectral tailoring. This 
is also required to obtain a long lived periodic structure in the ground state: 7s in 
our case critically depending on the temperature [22]. The AFC structure is built by 
optical pumping between the spin doublets. Although small (typically 2% for our field 
orientation [24]), the optical branching ratio is compensated by a long lifetime of the 
ground state population. This is suffiecient to efficiently pump the ions between the 
long-lived spin levels. The laser polarization is also applied along the [001] axis. In this 
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configuration the splitting is the same for the different crystallographic sites interacting 
with the laser. The bandwidth of the AFC is actually larger than A e because we use a 
moderate magnetic field. This could be an issue in practice because one would obtain 
the interference of two AFCs. In order to avoid this effect, we simply make the comb 
period 1/T coincide with the excited level splitting (A e = 2/T in our case). 
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Figure 1. Relevant level structure of the trivalent thulium under magnetic field. The 
optical pumping between the two Zeeman states is mainly allowed by a weak branching 
through the 3 H4(0) excited state. We also represent schematically the pulse sequence 
(see text for details). 



The spectral preparation is done by applying a sequence of pulse pairs (FWHM 
duration ~ 300ns separated by T = 1.5/is and then followed by a 100/xs dead time 
before the next pair). This train consists of 5000 pairs followed by a long waiting time 
(50ms) to ensure the complete decay of population from the 3 H 4 (0) excited state (lifetime 
~ 800/is). Each pair whose area is much smaller than it creates a small population 
difference [18] with a periodicity of 1/T in frequency space. After accumulation, the 
sequence gives rise to the atomic frequency comb. To first probe the frequency spectrum 
and to later generate a pulse to be stored, we use a very weak independent signal beam. 
Temporal shaping and frequency scanning is achieved by two independent acousto-optic 
modulators (AOM). The pumping beam has been also designed to be twice as large as 
the signal in the crystal. The AFC is then spatially uniform over the entire signal beam. 



2.2. Atomic Frequency Comb spectrum and storage efficiency 

By sweeping the signal AOM driving frequency, we can directly observe the transmission 
spectrum (see fig. [2^ and [2b). Even if the pumping power is very weak, we see that the 
AFC structure depends strongly on the preparation pulse intensity. For an increasing 
power from [2^ to [2b, the peak maximum absorption decreases due to power broadening 
of the hole-burning process (see insets). The sweep rate (320kHz) 2 being slightly too 



Efficient light storage in a crystal using an Atomic Frequency Comb 



5 




-2-1012 -2-1012 

Frequency (MHz) Frequency (MHz) 




-10 12 3-10123 

Time (|a,s) Time (\is) 



Figure 2. a) and c): AFC prepared for different pumping powers. The 660 kHz period 
is given by the delay between the preparation pulses T. The insets shows the central 
part of the optical depth spectrum. b) and d): A weak signal is sent matching the AFC 
bandwidth of a) and c) respectively. After being partially absorbed, it gives rise to a 
retrieval. The efficiency strongly depends of the shape of comb, which in turn depend 
on the power of the preparation pulses: the power is larger for the set c)-d) than for 
a)-c). 

high, we observe a distortion of the observed spectrum (see fig. [2b for example). This 
effect is not apparent in the insets representing only the central part of the spectrum 
where a slower sweep rate (120kHz f has been used. 

A weak pulse to be stored is then sent onto the medium. To be sure that the signal 
is uniformly covered by the comb, its duration is chosen to be slightly longer than the 
preparation pulses (FWHM ~ 450ns). When a signal is absorbed into the medium, it 
produces a retrieval. Curve [2b (resp. [2bl) corresponds toEK (resp. |2b). To be able to 
compare the incoming pulse and the retrieved signal, we simply bleach the medium by 
applying a strong single pumping beam (stronger than the preparation pulses). This 
procedure gives a reference for a fully transmitted incoming pulse. We can then estimate 
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the total efficiency: r] exp = 9.13 ± 0.10% (resp. r] exp = 4.37 ± 0.09%) for the comb [2b 
(resp. for[2ji). This corresponds to more than one order of magnitude improvement as 
compared to the previous proof of principle realization |18j . 

3. Discussion 

We now investigate the storage efficiency thanks to the independently measured AFC 
spectrum. To do so, we have developed an original model to analyze our data. The 
expected efficiency has already been calculated by assuming that the comb is composed 
of well-separated gaussian peaks [17]. This calculation was based on a temporal image 
by following the evolution of the atomic variables during the time sequence. Here we 
derive a more general formula valid for any peak shape. 

3.1. Expected efficiency 

We consider the propagation of a field £ in the spectral domain: 



The action of the medium is represented by the susceptibility x( a; ) is the wave vector). 
To account for the periodic structure of the AFC (period 1/T), we naturally decompose 
x(a>) in a Fourier series: 



The components corresponding to p < are dropped for causality reason. We here 
calculate the efficiency in the forward direction but our analysis could be extended to 
describe the backward retrieval [T7J. The field S(z,u) has the same periodic structure: 



Each term of the decomposition gives the amplitude of a retrieval centered at the time 
pT. The first term p = represents the absorption of the incoming signal. In the 
forward direction, the efficiency of the AFC protocol is then directly given by | a\ \ 2 . 
All the coefficients can be recursively determined by solving the propagation equation 
for a p (z). The first two components are given by: 

d 2 a (z) + k 2 (1 + c ) a (z) = 

d 2 z ai(z) + k 2 (1 + c ) at(z) + k 2 cia (z) = 

The efficiency in the forward direction is then given by 



d is the average absorption depth defined as d = — klm (cq) L where L is the length of the 
medium. Our formula is very general and just assumes that the susceptibility is periodic. 



d 2 £(z, u) + k 2 (l + x (oo)) £(z, u) = 



£(z, uS) = £(0, to) 2J a p (z) exp {—ImpuT) 

p>0 




(1) 
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It is sufficient to estimate the expected efficiency for a given frequency comb as presented 
in fig. [2j Our model doesn't make any assumption on the exact shape of the profile. This 
is the main advantage of our model. It also represents an alternative interpretation of 
the phenomenon in the spectral domain: the periodic filter is considered as a diffraction 
grating. A simple numerical Fourier expansion of the absorption spectrum (see insets 
of fig. [2^ and [2b) gives the different coefficients to determine rj. We calculate r] = 7.8% 
(resp. r\ = 3.5%) for fig. [2^ (resp. [2b). More generally for different preparation pulses 
intensities, we find a correct agreement between the experimental efficiencies and the 
ones deduced from the measured absorption spectrum (r] exp is represented by a red 
dashed line and i] by a black dashed line in fig. [3]). Nevertheless we systematically 
observe that the measured values are larger than the prediction. This may be due to 
inaccurate measurement of the absorption spectrum. As already mentioned, the lower 
sweep rate (120kHz) 2 used to obtain the central part of the spectrum may still be 
too high as compared to the width of the AFC peaks. This actes as a convolution 
effect, artificially broadening the peaks. More generally it is not fully obvious to 
accurately measure an optical thickness much larger than one. Any very far off resonance 
contribution of the laser would be transmitted through the medium. Here it could be due 
to the additionnal tapered-amplifier whose fluorescent background is extremely broad. 
For both effects, the real optical thickness would be underestimated and so would be 
the predicted efficiency. 

3.2. Current limitations of the pumping scheme 

It is now important to understand the current limitations. We will first estimate the 
maximum efficiency that could be obtain in our case and compare it to our measurements 
(see fig. [3]). The total efficiency strongly depends on the exact shape of the frequency 
comb. For a given initial optical thickness, one has to properly choose the peak width 
to optimize the retrieval. We now assume that the comb is composed of lorentzian 
(half width at half maximum V) periodically distributed peaks. The lorentzian shape 
allows simple analytical calculation and gives a satisfying agreement with the measured 
spectrum. The efficiency depends only on the maximum optical depth d and the comb 
finesse F = n/ (TT) p2]. 



This result can be compared to eq. 1 of [THj. The structure and the term 
interpretation are similar. Our formula includes by definition a potential overlap of 
the lorentzian peaks when the finesse is low. This has been alternatively interpreted by 
de Riedmatten et al. as a uniform absorbing background and fitted independently [18J. 
In our case, we can directly fit the AFC spectrum with eq. [2j Our result can be also 
simplified when the finesse is high. In that case, the best efficiency of the protocol is 
expected. For a given initial optical thickness d, the finesse should be properly chosen 
to obtain a optimized efficiency. For F opt — tt (1 + d/A), we have 



r? = ci 2 tanh 2 (tt/2F) e -^H-/m e -^/F 



(2) 



Vopt = 4e~ 2 d 2 / (4 + df 



(3) 
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The efficiency is limited to 54% in the forward direction. This value is a common feature 
of the protocols involving a dipole rephasing [15]. 

One can now wonder if the measurements correspond to the optimum for a given 
optical depth. As we have seen, the shape of the comb is essentially controlled by varying 
the intensity of the preparation pulses. At the same time this procedure is changing the 
maximum absorption depth d and the width of the peaks (the finesse F). For different 
pumping power, we can directly measure the maximum absorption depth d and then 
estimate the corresponding r] opt (eq. [3]). 




Figure 3. Color online: Measured efficiency i] exp as a function of the maximum 
optical depth (red dashed line used to guide the eye). We have also calculated the 
expected efficiency r\ (eq. [TJ from the absorption spectrum (black dashed line used to 
guide the eye). Then assuming that the comb is composed of lorentzian peaks (eq. 
[2j, we have finally plotted the optimum efficiency rj opt (solid black line) and fitted 
the experimental spectrum to obtain the finesse. Inset: Measured finesse obtained by 
fitting the absorption spectrum with a comb of lorentzian peaks (red dashed line used 
to guide the eye) and comparison to the optimum F opt (solid black line). The two 
cases presented in fig. [2] are circled. 



The comparison between the measured efficiencies r] exp (red dashed line in figj3]) 
and the expected optimum for a given optical depth i] opt (solid black line in figjSJ) allows 
us to evaluate the quality of the preparation. There is clearly a good matching of the 
efficiency curves at low optical depth (between and 3). In that case, the efficiency is 
not critically depending on the finesse (eq. EJ) [lZ]. In practice any finesse between tc 
and 2tt could be considered as optimal. This is verified experimentally and explains the 
efficiency curve matching even if the finesse is only roughly following F opt (see inset in 
figJ3|). On the contrary, we observe an efficiency breakdown at larger optical thickness 
(d ~ 3 in figj3]) where r\ is now critically depending on the finesse (eq. |2J). The measured 
finesse is indeed significantly deviating from F opt (see inset in figHS). This is sufficient to 
explain the efficiency breakdown in this region. 
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This clearly shows that the preparation sequence may not be appropriate at large 
optical thickness. Even if we should be able to benefit from the initial optical depth of 
our crystal (d ~ 4 — 5), we cannot properly prepare the AFC to obtain the maximum 
efficiency In practice the external control of the preparation is given by the power of 
the pumping pulses. This single degree of freedom cannot independently master the two 
parameters of the comb (height d and finesse F). Our comparison shows that a simple 
pumping strategy using a train of pulse pairs is relatively appropriate at low optical 
thickness as soon as the power is properly adjusted. On the contrary, at large optical 
thickness, the efficiency is limited. 

The pumping intensity primarily affects the width of hole-burning process because 
of the power broadening. This effect is intimately related to the optical pumping 
dynamics. To break this limitation, it is possible to apply a coherent population transfer 
scheme [251 [261 EZj - This approach has been successfully implemented to obtain and 
isolate narrow feature in different systems including Tm 3+ :YAG [28l l26l [27] . In other 
words, by enriching the preparation sequence with numerous pulses, one would obtain 
another degree of freedom to control accurately the comb shape. The propagation of 
sophisticated pulse sequence in highly absorbing materials may nevertheless limit the 
capability of the coherent population transfer. In that case, one can also use a transverse 
illumination of the slab [29]. Combining this two techniques, it should be possible to 
independently control the efficiency and the spectral properties of the AFC. 

4. Storage at the single photon level 

Even if the efficiency is a critical parameter for a quantum storage protocol, a noise 
estimation should also be done. In order to investigate this level in the very weak field 
condition, we reproduce the experiment with less than one photon per pulse. 

We use the same preparation procedure that have been optimized previously. The 
weak incoming pulse to be stored is strongly attenuated by using neutral density filters. 
After transmission through the crystal, the signal is now fiber-coupled to a single photon 
counting module (Perkin-Elmer SPCM AQR 14 FC). To increase the measurement 
statistics, we repeat 1000 signal pulses after each preparation procedure (5000 pairs 
followed by a 50ms waiting time). The total duration of the sequence is now 329ms 
corresponding to 3039 signal pulses per second. The arrival of the photon is counted 
in a 300ns time gate and accumulated for 5.51 seconds. The corresponding histogram 
is presented in fig. H] for different position of the time gate covering the arrival of the 
transmitted signal and the AFC echo. 

We have independently measured the collection efficiency from the crystal to the 
single photon counting modules. Knowing the quantum efficiency of the detector and 
assuming that the AFC efficiency is the same as previously measured, we estimate the 
incoming number of photons per pulse in the crystal to be ~ 0.5 photon/pulse. 

For comparison, the measurement has been also done at higher power (solid curve 
in fig. 0J. We observe the same kind of trace as previously described (see fig. [2] for 



Efficient light storage in a crystal using an Atomic Frequency Comb 



10 



120 




-2-10123 



Time (|is) 

Figure 4. Color online: Histogram of the accumulated counts for 5.51s within a 
300ns time gate. The incoming pulse contains ~ 0.5 photon/pulse. We have recorded 
independently the signal at higher power using an standard photon diode (solid red 
curve). The three observed pulses are respectively the transmitted power, the retrieved 
signal (first echo) and a second echo (see text for details). 

example). Nevertheless we suspect the efficiency to be higher because we also observe a 
second echo. This can be due to the slightly different experimental conditions (magnetic 
field). We haven't measured the efficiency in that case but this effect is a potential 
route to improve the retrieval and is currently under investigation. The photon counts 
(histogram in fig. |4|) follow the same curve within the statistical fluctuations. We 
additionally observe a detection background of ~ 8 counts. By comparing this noise 
level to the maximum of the retrieved signal, we obtain a signal to noise ratio of 11. The 
intrinsic dark counts of the detector should be of 0.31 count with this cycling, gating 
rate and accumulation time. We have isolated this detection background. It is due to 
a leak of the preparation beam. Even if we use a first AOM to temporally shape the 
preparation pulse and a second one right before the detector, the total isolation is not 
sufficient to completely remove the residual light. In the near future, we can also use 
a mechanical chopper since the waiting period is long (50ms). We haven't finally seen 
any spurious light due to the medium itself (fluorescence). The spontaneous emission is 
strongly reduced after a long waiting time and is also emitted in all directions. We then 
believe that the ultimate limitation of such a measurement are the dark counts of the 
detectors. In other words, the AFC storage should not add any noise to the incoming 
photons. This is a remarkable figure of merit that make the protocol promising for the 
storage of non-classical field. 

5. Conclusion 

To conclude, we have successfully applied the AFC protocol to a Tm 3+ :YAG crystal 
whose properties are of particular interest in the prospect of quantum repeaters. We have 
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observed efficiencies up to 9.13 ±0.10%. Our measurements correspond to more than an 
order of magnitude improvement as compared to previous realizations in Nd 3+ :YVC>4 
[18J. Since the material properties gives the ultimate width of the spectral selection, 
they strongly influence the efficiency. This comparison shows that further improvements 
will demand material developments. This is especially true for the future backward 
configuration since Tm 3+ :YAG and Nd 3+ :YVC>4 have only two ground states. This 
configuration does not allow, at the time, optical pumping in an auxiliary level and the 
application of control pulses [TT] . 

We have developed an original spectral analysis. Our model directly predictd 
the observed efficiency for a given initial absorption spectrum that we independently 
measure. For a definite bandwidth and a fixed delay, the efficiency is both limited by 
the available optical thickness and the preparation method. In that case, due to power 
broadening effect, the width and the height of the comb peaks are correlated. As a 
consequence, the efficiency and the spectral feature of the protocol (delay and multi- 
mode capacity) are not independent. These characteristics are critical for quantum 
repeater applications. A global estimation of the ultimate performances will demand 
the development of original preparation methods. 

We have finally shown that even for very weak pulse ~ 0.5 photon/pulse, the storage 
itself should not add any noise to the incoming field. The efficiency and the noise feature 
render then the AFC particularly attractive as a quantum memory protocol. 

References 

[1] H. Briegel, W. Duer, J.I. Cirac, and P. Zoller, PRL, 81, 5932 (1998). 

[2] B. Julsgaard, J. Sherson, J. Fiurasck, J. Cirac, and E. Polzik, Nature 432, 482 (2004). 

[3] M. Eisaman, A. Andre, F. Massou, M. Fleischhauer, A. Zibrov, and M. Lukin, Nature 438, 837 

(2005) . 

[4] T. Chaneliere et at, Nature 438, 833 (2005). 

[5] K. S. Choi, H. Deng, J. Laurat, and H. J. Kimble, Nature 452, 67 (2008). 

[6] Y.-A. Chen, S. Chen, Z.-S. Yuan, B. Zhao, C.-S. Chuu, J. Schmiedmayer, and J.-W. Pan, Nat. 
Phys. 4, 103 (2008). 

[7] M. Atature, J. Dreiser, A. Badolato, and A. Imamoglu, Nat Phys 3, 101, (2007). 
[8] L. Childress et at, Science 314, 281 (2006). 

[9] E. Fraval, M. J. Sellars, and J. J. Longdell, Phys. Rev. Lett. 95, 030506 (2005). 
[10] J. J. Longdell, E. Fraval, M. J. Sellars, and N. B. Manson, Phys. Rev. Lett. 95, 063601 (2005). 
[11] A. Louchet et at, Phys. Rev. B 77, 195110 (2008). 
[12] M. Nilsson and S. Kroll, Opt. Commun. 247, 393 (2005). 

[13] A. L. Alexander, J. J. Longdell, M. J. Sellars and N. B. Manson, Phys. Rev. Lett. 96, 043602 

(2006) . 

[14] B. Kraus, W. Tittel, N. Gisin , M. Nilsson, S. Kroll, J.I. Cirac, Phys. Rev. A 73 020302(R) (2006). 
[15] G. Hetet, J. J. Longdell, A. L. Alexander, P. K. Lam M. J. Sellars, Phys. Rev. Lett. 100, 023601 
(2008). 

[16] W. Tittel, M. Afzelius , R. L. Cone , T. Chaneliere , S. Kroll , S. A. Moiseev , M. Sellars, 
larXiv:0810.0172l 

[17] M. Afzelius, C. Simon, H. de Riedmatten, and N. Gisin. larXi"v:0805.4164l 



Efficient light storage in a crystal using an Atomic Frequency Comb 



12 



[18] Hugues de Riedmatten, Mikael Afzelius, Matthias U. Staudt, Christoph Simon and Nicolas Gisin, 
Nature, 456, 773 (2008). 

[19] J. Nunn, K. Rcim, K. C. Lec, V. O. Lorenz, B. J. Sussman, I. A. Walmsley, and D. Jaksch, Phys. 

Rev. Lett. 101, 260502 (2008). 
[20] Christoph Simon et al, Phys. Rev. Lett. 98, 190503 (2007) 
[21] R. M. MacFarlane, Opt. Lett. 18, 1958 (1993). 

[22] N. Ohlsson, M. Nilsson, S. Kroll, and R.Krishna Mohan, Opt. Lett. 28, 450 (2003). 

[23] B. Lauritzen, S. R. Hastings-Simon, H. de Riedmatten, M. Afzelius, and N. Gisin, Phys. Rev. A 

78, 043402 (2008). 
[24] F. de Seze et al, Phys. Rev. B 73, 085112 (2006). 
[25] I. Roos and K. M0lmer, Phys. Rev. A 69, 022321 (2004). 

[26] F. de Seze, F. Dahes, V. Crozatier, I. Lorgere, F. Bretenaker and J.-L. Le Gouet, EPJD 33, 
PP 343-355 (2005). 

[27] A. L. Alexander, R. Lauro, A. Louchet, T. Chancliere, and J. L. Le Gouet, Phys. Rev. B 78, 
144407 (2008). 

[28] L. Rippc, M. Nilsson, S. Kroll, R. Klieber, and D. Sutcr, Phys. Rev. A 71, 062328 (2005). 
[29] M. Tschanz, A. Rcbane, D. Reiss, and U. P. Wild, Mol. Cryst. Liq. Cryst. Sci. Tcchnol., Sect. A 
283 (1996) and reference therein. 



